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I. INTRODUCTION 



From its inception, lattice gauge theory was hampered by the absence of an adequate 
fermion formulation. It was thought for many years that there is no discretization scheme 
that displays explicit chiral symmetry without either sacrificing locality or introducing non- 
physical doublers in the fermion spectrum. This state of affairs was summarized in the 
famous no-go theorem of Nielsen and Ninomiya 0. Recent years have seen major progress 
in the construction of lattice regularizations of fermions with good chiral properties J2|. 
Domain wall and overlap 0] formulations preserve chiral symmetry exactly, allowing the 
construction of vector gauge theories on the lattice. An idea common to these schemes is 
that of employing a large number of fermion flavors, one of which survives in the continuum 
limit as a massless fermion supporting exact chiral symmetry. 

In the domain wall formulation, the expanded flavor space may be seen as an extra 
dimension with a defect in a background field. One chirality of the Dirac spinor is expo- 
nentially localized along the defect. If the extra dimension is periodic and of finite extent 
L 5 , the other chirality will be localized along an unavoidable second defect, while all the 
other Dirac fermions, L5 — 1 in number, stay heavy. This is also the scheme in Shamir's 
surface-fermion variant |||6|], which we study in this paper; here the chiral modes are local- 
ized on the surfaces of a five- dimensional slab. In either case, the limit L 5 — > 00 requires 
that the contribution of the heavy flavors be subtracted by introducing L 5 flavors of 

pseudofermion fields, often called Pauli-Villars (PV) fields. These are bosons that have the 
same index structure as the fermions. 

Domain wall fermions are equivalent to overlap fermions Q in the limit that the size and 
lattice spacing of the fifth dimension are taken to infinity and zero, respectively [0,0. The 
bulk degrees of freedom remain massive. They decouple as the four-dimensional lattice is 
removed, leaving a four- dimensional Dirac operator describing the chiral surface modes. 

In this paper we study domain wall fermions in the limit of strong gauge coupling. We 



adapt the classic method of Kawamoto and Smit [I0| to derive an effective action for mesonic 
degrees of freedom.^] A novel ingredient of the domain wall formalism is the PV bosons. 
We combine the fermions and pseudofermions into a supervector, and similarly group the 
effective degrees of freedom into a supermatrix. Then the formalism of supergroups allows 
us easily to extend the treatment of the fermions to an analysis of the complete theory. 
The effective action contains all the meson-like states that can be constructed from pairs of 
fermions, from pairs of bosons, and from boson-fermion pairs. 

The interesting question is, does this theory really behave as one would expect of a theory 
of exactly chiral quarks? To address this, we take the number of colors N c to infinity and look 
for saddle points of the effective action. Following Kawamoto and Smit, we begin by setting 
the Wilson parameter r to zero and classifying the symmetries of the theory; we display 
explicitly the Goldstone bosons that are the fluctuations around the saddle point. The 
spectrum of Goldstone bosons is characteristic of a theory with complete fermion doubling, 



1 See also [O]. This method has been elaborated and applied many times. See for instance 
lEMLj. 
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which is not surprising since it is the term proportional to r in the fermion action that 
breaks naive doubling in the first place. When we proceed, however, to perturb the saddle 
point with the r terms, we find that the doubling persists and that the only effect of the 
perturbation is identical to that of adding a mass term (a la Shamir to the domain wall 
action. We conclude that in the strong coupling limit, domain wall fermions are doubled 
and massive, with all that that implies for the mesonic spectrum. 

We thus confirm the conclusion reached by two of us |TB[ through study of the Hamilto- 
nian of domain wall fermions at strong coupling. Unfortunately, our analysis here is more 
involved and perhaps less conclusive. Study of the Hamiltonian via Rayleigh-Schrodinger 
perturbation theory is straightforward and leads unambiguously to an effective Hamiltonian 
for mesonic degrees of freedom. The symmetries of this effective theory are readily appar- 
ent. In the Euclidean formulation some additional approximation is necessary, such as a 
hopping parameter expansion, mean field theory, or the large- N c limit .0 Even at large N c , 
one might argue that there exists another saddle point with different symmetry properties. 
Nonetheless, the Euclidean study complements the Hamiltonian results in that it does not 
require a time-continuum limit and thus does not introduce a new question regarding the 
order of limits.^ 

In Sec. [TT] we derive the strong-coupling effective action. For pedagogic reasons, we 
begin by ignoring the PV bosons. Following Kawamoto and Smit |fT0| , we integrate out 



the gauge field in this limit and rewrite the resulting path integral in terms of a bosonic 
matrix field. That this path integral generates the same Green functions as the original 
Grassmann integral is due to an equality between generating functions that was proven in 
|T0| (following [[HJ)- This path integral contains the effective action. While its derivation 
does not depend on a large- iV c limit, we give its explicit form for this case only. In the 
last part of Sec. |TJ we repeat the derivation with the bosons included. This requires a 
generalization to supermatrices of the identity relating generating functions, which we prove 
in Appendix [B| 

Since the chiral symmetry of the domain wall action arises dynamically, it is difficult 
to analyze the symmetry of the effective action directly. The symmetry analysis cannot 
be divorced from the dynamics. Section JT| is devoted to deriving the symmetry breaking 
pattern of the theory and the spectrum of Goldstone bosons. We can only make analytic 
progress in the N c — > oo limit. Here the task at hand is to find saddle points of the effective 
action. Again we start with an analysis of the theory without the PV bosons. We begin 
by setting r to zero which makes finding a saddle point straightforward. In this limit the 
symmetry of the effective action is apparent; upon introducing an ansatz for the saddle point 
we display the pattern of spontaneous symmetry breaking and the concomitant Goldstone 



In the Hamiltonian formalism, a strong coupling expansion combined with a large- N c expansion 



has been used in [|TJ],[lj| . 



3 It is curious but true that one may consistently ignore the PV bosons in the Hamiltonian 
calculation to lowest order in 1/g 2 [16|, but not in the Euclidean calculation. To our knowledge, 



there is no established correspondence between the perturbation series in the two formalisms. 
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bosons. Then we turn on the Wilson term in the action and we compute its effect on the 
saddle point. It turns out that the degeneracy of the r = saddle point is broken by the r 
term, but only in the way that a mass term breaks it, that is, by breaking chiral symmetry 
while leaving the doubling. We show explicitly that inclusion of the Shamir mass term does 
the same. 

The extension of the calculation to include the PV bosons is straightforward. The saddle 
point in the enlarged field space is the obvious extension of the saddle point found above. 
Perturbing the saddle point with the Wilson term of the bosonic fields gives a mass matrix 
for the pseudo-Goldstone bosons that is identical in form to that due to the massive fermions. 

We conclude with some discussion and comparison to other work. 



II. THE EFFECTIVE ACTION AT STRONG COUPLING AND LARGE N c 

A. Without Pauli-Villars bosons 

We define domain wall fermions on a five-dimensional hypercubic lattice. The fifth 
dimension has finite extent L5; for generality we assign a lattice spacing a to the four 
dimensions of Euclidean space-time and a separate lattice spacing 05 to the fifth dimension. 
We shall denote by x = (x\, X2, £3, £4) a coordinate vector in the four dimensional space, 
with Greek indices /x, u, . . .; the fifth coordinate is s. The fermion field ip^ * s a Grassmann 
field carrying color index a = 1, . . . , N c and Dirac-flavor index a = 1, . . . , 4Nf. The gauge 
field U Xfl is an SU(N C ) matrix field that resides on the four- dimensional links only and is 
independent of s. 

The action is a sum of pure-gauge and Dirac actions, 

S(U,$,1>) = Su(U) + S F $,1>,U). 
The gauge action is as usual a sum over plaquettes, 



S, 



We write the fermion action as 



(2.1) 
(2.2) 



witl 



Sf — S4 + £5 + S m , 



XS[l 



aa 
xs 



a/3 



Ls-1 r 



Sm, 



aa 
xs 



1+75 



x s=l L 

(m-4-p)£OA 



a/3 



+ t 



aa 

x,s+l 



r -7/* 

2 / a /3 
1-75 



aft 



aa 

xs 1 



(2.3) 

(2.4) 

(2.5) 
(2.6) 



4 Our Dirac matrices 7 M are hermitian and satisfy {7^,7^} = 25^ u . We define = 74 7 M so that 
04 = 1 and 01,2,3 are anti-hermitian. 75 is hermitian. (See Appendix |A|.) 



and p = a/a*,. Only S4 contains the gauge field. 

In the strong-coupling limit we drop Sjj. The effective action will then be written in 
terms of the meson field, 



M a £(x 



N, 



(2.7) 



c a=l 



This is a color singlet, local in x but nonlocal in s, and hence we write it as a (4NfL 5 ) x 
(4NfL 5 ) matrix. (This is in line with the view that s represents an internal flavor index. 
Eventually we will isolate the projection of M. onto the subspace spanned by the chiral 
surface modes.) 

In the strong-coupling limit the partition function is given by 



V^VifjVU expS F (ip,ip,U) 



(2.8) 



The integral over the gauge field can be performed immediately by noting that it is simply 
a product of decoupled link integrals, 



VU expS 4 (ifj,ij,U) 
YlJdU exp Tr [A^x)U + U ] A^x% 



where we have defined 



a/3 



1>: 



b/3 



A b ;( x ) 



act 

x+ji,s 



a/3 



ris 



(2.9) 

(2.10) 
(2.11) 



The integral in Eq. (|2.9|) can be carried out explicitly for any value of iV c PU| , and the result 



expressed in terms of M.(x). We are interested in the limit N c — > 00 [21 , |22| . Defining 

X^x) = -4M(x) T -^^M{x + jl) 7 -^, (2.12) 



the result is 



S 1 (M) = N c J2trF(X^x)), 



XfJ, 



where 



F(X) = 1 - y/l-X + log 



1 + VT^X 



(2.13) 



(2.14) 



The trace ("tr") in Eq. ( |2.13| ) denotes a trace over Dirac-flavor and s indices, not to be 
confused with the color trace ("Tr") appearing in Eq. ( |2.2| ). We plot F(X) in Fig. [I], and 
note that one may use a straight line F(X) ~ A/4 as a first approximation. 



5 



The remainder of Sf couples fields along the fifth axis for fixed x. We can rewrite this 
in terms of M. as well, viz., 



s ± = s 5 + s 



M 



x) 



where we denote by D 1 - the fifth-dimension Wilson-Dirac operator, 

75* 



Dts> = P 



1 + ,/l 

S ,s'-l + 



S , S '+1 



+ (m - 4 - p)8 a 



(2.15) 



(2.16) 



2 / ' \ 2 

with implicit limits on s and s' as indicated in Eq. fl2.5|) . The complete partition function 
is then 

(2.17) 



Z(J) = I V^Vij; exp [S^M) + S ± (M) + Sj(J,M) 
where we have added an external current coupled to mesonic bilinears according to 



Sj = £4W(*)« = -N c J2trJ(x)M( 



x 



(2.18) 



so that functional differentiation with respect to J generates Green functions of the fermion 
bilinear field A4(x). 

As shown in JlO| , the same Green functions are generated by the purely bosonic partition 
function 



Z(J) = / VM(x) exp 



-iV c £tr logM(x) 



exp 



5i(M) + S L {M) + Sj(J,M) . (2.19) 



Here M(x) G U (ANfL 5 ) is a bosonic, unitary matrix field and T>M(x) indicates the invariant 
Haar measure at each 4d site x. (We derive this result in Appendix [FJ, and generalize it to 
include Pauli-Villars fields.) Collecting everything, we have 



S C «{M) 



with 



S eS (M) = N c 



Z(J) = J VM(x)e 
tr F(-M{x) (r + ^)M(x + p)(r - ^] 

X,fJ, 

-J2 tT lo S M ( x ) ~J2 tr J(x)M(x) 



[2.20) 



x) 



(2.21) 



X X 

We set J = henceforth. 

The fifth-dimension Wilson-Dirac operator D x is not hermitian. It will be convenient to 
work with the hermitian operator h = 74-D -1 , so we substitute M — > M74, giving 



S eS (M) = N c 



tr F(-M(x) lA {r + 7 M )M(x + fl)j 4 (r - 7m )) - E tr hM ( x ) 

X,fl X 



£tr logM( 



x 



(2.22) 



The measure and tr log M are unchanged since 74 is a special unitary matrix, h is precisely 



the single-site Hamiltonian studied (in the — > limit) in | 16| . We review its spectrum 
and eigenf unctions in Appendix 



B. Addition of Pauli-Villars bosons 



We have so far neglected the Pauli-Villars fields. Introduction of supermatrix notation 
will enable us easily to extend our calculation to include their effects in the effective action. 

Following Vranas ||, we introduce scalar fields 0"" and their conjugates, with the same 
indices as the fermions. Their action Spv{4>,4>) is identical in form to the fermion action 
Sf [see Eqs. (|2.3|) - (|2.6|) 1. The only difference is the imposition of antiperiodic boundary 
conditions instead of the free surfaces at s — 1 and s = L5, which we write as a modification 
of S 5 , 



u r 



pEE 



X s=l . 



laa 
xs 



1+75 



a/3 



75 



a/3 



xs 



where 



(2.23) 



Laa 



We group ip and 4> into a superfield, 



xs 



( J,aa 

TIS 



V 



Laa 
xs 



(2.24) 



(2.25) 



The generalized meson field is a supermatrix containing Grassmann elements as well as 
ordinary numbers, 



iVe 



a/3 



N c ( J.aa i a(3 i aaja/3 \ 
VxsVxs' VxsVxs' 



N, 



c a=l 



N, 



c a=l 



(2.26) 



It is a {ANf2L^} x {ANf2L^) matrix that contains the old meson matrix ( |2.7| ) in the upper 
left-hand corner. 

In the strong coupling limit the partition function is now 



VipVi)V<f>V(f>VU exp [S F (iJ),il>, U) + S PV {<p, 0, U) 
We integrate out the gauge fields as we did in Eq. ( |2.9| ), 

e s^M) = JvU expfS^, ^ U) + SS, 0, U)] 
= l[[dU exp Tr [A^x)U + E7%(a;)], 

where we have defined, in analogy with Eqs. fl2.10| ) and (2.11), 

'r + 7/ 



(2.27) 



A b ;(x) 



E* 

s 

E* 



aa 
xs 



a/3 



aa 

x+fi,s 



a/3 



(2.28) 

(2.29) 
(2.30) 
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Performing the integral (|2.28|) in the limit N c — > oo we obtain [cf. Eq. fl2.13|) 1 

S 1 (M)=N c J2^TF(\ fl (x)), 



(2.31) 



where Str denotes the supertrace (see Appendix |B|). F and A are the same as before, with 
A defined in terms of the new supermatrix field Ai. The equivalent of Eq. fl2.15| ) is 



S 1 = -N C Y, Str D ± M( 
where the new D L is defined by 

D 1 - 



x 



(2.32) 



D F 
D 



(2.33) 



PV 



Here Dp is the fermionic fifth-dimension Wilson-Dirac operator given in Eq. (|2.16| ); the 
Wilson-Dirac operator Dp V for the bosons is derived from Eq. ( p.23| ) and differs from Dp 
only on the boundaries. Thus we arrive again at Eq. fl2.17|) , 



P^D^exp S 1 (M) + S ± (M) 



(2.34) 



mutatis mutandis. 

We show in Appendix [B| that the integral over the superfield \1/ can be replaced by an 
integral over a supermatrix field M, 



VM(x) exp 



-AT c ^Str logM( 



X) 



exp 



S 1 (M) + S ± {M) 



(2.35) 



In this equation M(x) is an element of the supergroup U(4:NfL 5 \4:NfL 5 ) and T>M(x) is the 
corresponding invariant Haar measure at each 4d site x. Since D x is not hermitian, we again 
substitute M — > M74 and define the hermitian operator h = 74-D- 1 -. The result is 



S eS (M) = N c 



£Str F(-M(x) 74 (r + l„)M{x + //) 74 (r - 7/ j) - E Str ^f 1 ) 



-^Str logM(x) 

X 

This is the effective action for all the meson-like degrees of freedom. 



(2.36) 



III. CHIRAL SYMMETRY BREAKING PATTERN AND GOLDSTONE BOSONS 

The large- N c limit allows us to evaluate the partition function by finding saddle points of 
the effective actions ( |2.22| ) and ( |2.36| ). We shall begin by setting r = 0, whereupon Wilson 



fermions reduce to naive fermions, and then we will see how things change as the Wilson 
term is restored. Again, we leave the inclusion of PV bosons to the end. 
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A. The r = case; no PV bosons 



For r = the effective action fl2.22| ) becomes 
S eS (M) = N c 



x 



(3.1) 



^trF(M(x)a M M(x + /t)a M ) - ^tvhM{x) -£}tr logM(; 

_X,fi X X 

We can eliminate the a matrices by performing the following unitary transformation (spin 

(3.2) 



diagonalization |T7|j23]| ) 

M(x) - {a^i^T^a^Mix^ri^ria^. 
We obtain the simplified expression, 



S eS (M) = N c 



tr F(M(x)M(x + /})) - J2 tr hM(x) -J^tr log M(x) 



(3.3) 



If we set h = 0, the theory described by the effective action (|3.3|) is invariant under the 
transformation 



M(x) 



{UM{x)V^ for x even 
VM(x)W for x odd 



(3.4) 



where U, V G U(4NfL^). This is a global left-right symmetry typical of non-linear sigma 
models. Inclusion of h restricts the transformation by the condition Uh = hV. Within 
a degenerate eigensubspace of h, this means that U — V. Within the nullspace of h, the 
restriction vanishes and the symmetry under U ^ V rotations returns. In the limit L 5 — > oo, 
h possesses 4Nf zero modes, which are just the chiral surface modes for the zero-dimensional 
domain wall. The symmetry realized among these modes is U{ANf) x U(ANf). 

To translate this back into the original fermionic coordinates we must note that we have 
redefined M(x) through left multiplication by a factor of 74 so that the standard chiral 
U(Nf) x U(Nf) transformation is 



(3.5) 
(3.6) 



where the A a are U(Nf) flavor generators. This operation does not leave h invariant. The 
zero eigenvectors of h, however, are eigenvectors of 75 (see Appendix ^) and thus the 
nullspace of h is invariant under this U(Nf) x U(Nf) group. This is a subgroup of the 
U(4Nf) x U(4Nf) symmetry group of the nullspace, in fact, of its U{ANf) subgroup defined 
by U = V. 

The group U(4Nf) x U(ANf) is the symmetry of fully doubled, naive fermions. Evidently 
for r = the domain wall formalism yields ordinary doubling. 

These symmetry considerations guide us in choosing an ansatz for the saddle point of 
S c s. Let us consider the eigenstate basis of the single-site Hamiltonian h, 



h\e) = e\e). 



(3.7) 
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We suppose that, at the saddle point, M(x) is translation-invariant and diagonal in the h 
basis, 

M(x) = M = ^2\e)M e {e\. (3.8) 

The effective action (|3.3|) now reads 

S eS (M) = N C VY, [d F ( M D ~ l ogM e - eM e ] , (3.9) 

where d = 4 is the number of dimensions. Using the formula ( |2.14j ) for F, we differentiate 
with respect to M t to find the roots 



e(d - 1) ± Je 2 (d - l) 2 + (2d - l)(d 2 + e 2 ) 

Mf = - V , 2 2 -■ (3.10) 

d z + e 2 

In the subspace of zero modes of h, where e = 0, this reduces to 

M ± = ± y —^—- (3-11) 

We plot M ± (e) in Fig. ^|, and note that \M ± \ < 1 for any e. To maximize the action at the 
saddle point, we choose 

M,= \S f r£> ° (3.12) 
\ M e for e < 0. v 7 

This choice satisfies 

M_ e = -M e . (3.13) 

The solution ( |3.10| )-( |3.12| ) does not break any U — V symmetries [see Eq. (|3.4|) 1 within 
degenerate eigensubspaces of h, but it does break the U ^ V symmetries in the space of zero 
modes according to U(ANf) x U(4Nf) —>■ U(ANf). We may use a U ^ V transformation to 
set Mq = Mq without loss of generality, and so M = MqI within the space of zero modes. 

In order to display the Goldstone bosons, we allow M(x) to fluctuate about the saddle 
point M. We write 

and 

S(x) = e iH{x \ (3.15) 

We take the Goldstone field H(x) to be non-zero only within the zero modes of h, and thus 
it belongs to the algebra U(ANf); it commutes with M. Inserting Eqs. ( |3.14j ) and ( 3.15 ) 
into the action ( p.3|) , we expand to second order in H(x) to obtain[| 



5 The term linear in H that comes of expanding F cancels against the logarithmic term in Eq. 
In d = 4 the coefficient in square brackets in Eq. ( 3.16| ) is equal to 31/192. 
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S eS (M) ~ S cS (M) - i \F'(M 2 )M 2 + F"(M 2 )M 4 ] £ tr [H(x) - H(x + A)] 2 . (3.16) 



X/L4 



This is a quadratic action representing massless Goldstone bosons in the adjoint represen- 
tation of U(4N f ). 

If the field H is taken to represent fluctuations outside the U(ANf) subgroup, it acquires 
a mass of the order of m — 4, indicating that these bosons are not Goldstone bosons. We 
may sum up the results so far by stating that domain wall fermions with r = possess the 
symmetry of naive fermions and yield Goldstone bosons according to the simplest scheme 
of spontaneous symmetry breaking. 



B. The r / case; no PV bosons 



Turning on the Wilson term in the action will break the U(4Nf) x U(ANf) symmetry. 
In principle, we should redo the above analysis and find a new saddle point of the complete 
action. To simplify the calculation, however, we will stay within the manifold of saddle 
points of the r = action, and determine the effective action that selects the point of lowest 
energy in this manifold. The symmetry of this effective action will be the symmetry left 
unbroken by the r terms. 

We begin again with the action (|2.22|) , this time keeping r ^ 0. After spin diagonalization 
(p.2|), the strong-coupling effective action takes the form [cf. Eq. (|3.3|) 1 



S eS (M) = N c 



^ trF(M(x)[l - (-l) Xk n k ]M(x + k)[l - (-l) Xfc r 7fc ] 

x 

:=1,2,3 

f ^]trF(M(x)[l + s 3 (x)n 4 ]M(x + 4)[1 - s 3 (x)rj A 



- ^ tr hM(x) -J2 tr lo S M ( 



x) 



S^°(M) + AS eS (M). 



(3.17) 
(3.18) 



where S3 (or) 
which gives us 



AS eS (M) 



N c 



-l^i+^+a^ For the sake of simplicity we shall approximate F(X) ~ A/4, 
[-r(-l) Xh tr7 fc {M(x), M(x + k)} + r 2 tr M(x)-f k M(x + k)-f k 

x 

fc=l,2,3 

+ H [-rs 3 (x)tr^[M(x),M(x + 4)] - r 2 ti M(x)^ A M(x + 4)74 ). (3.19) 



We stay with the ansatz M(x) = M as given by Eqs. ( |3.8| )-( |3~T^ ), translation-invariant 
and diagonal in the eigenbasis of h. We leave free the orientation of M within the space of 
zero modes. The resulting effective action is 
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rr> 2 



AS eS (M) = N c - 



fc=l,2,3 



(3.20) 



We collect in Appendix [A| some formulas connected with diagonalization of h (see [TO). 
We work in the limit 05 — > 0, where /i becomes 

/i = 7475^ + (m-4) 74 . (3.21) 

In the limit L = L5CI5 — > 00, the 4Nf surface states of h become exact zero modes. We 
denote them (for each flavor) by (OV), where 77 = ± is the sign of the energy before taking 
L — > 00, and a is a spin label. We label states in the continuum spectrum of h as jeer); 
they satisfy h\ea) = e|ecr), with |e| > m — 4. It is straightforward to evaluate the traces in 
Eq. ( |3.20[ ) in the basis |0 ± cr) , |eer). Using the key formula ( |A16| ), we find that the first term 



in Eq. ( |3.20| ) is independent of the orientation of M in the zero modes (OV); the second 
term yields 



AS eS (M) = const. - N c — ^(0V|7 4 M7 4 |0 ? ' a'){tf> a'|M|0V). (3.22) 

^ X 

Only the second matrix element in Eq. ( |3.22| ) depends on the orientation of M in the zero 
modes. As we show in Appendix [A], the coefficient matrix ^alj^M'-f^O 11 a') takes the form 
Jpz in the indicated basis for the zero modes, where J is plotted in Fig. |3|. 

Now we change basis [16| in the space of zero modes from |0 ± cr) to \0 L,R a) = ^(|0 +cr ) =t 



|CT(t)). The new basis states are localized on either boundary and hence represent the true 
chiral modes. The change of basis converts ps to p\ (which is 74 in the chiral basis for the 
Dirac matrices). Hence the effective action is 



r 2 



ASeff(M) = const. - N c —jJ2^PiM. (3.23) 

^ X 

Alternatively, one may remove the p\ from Eq. (|3.23|) by means of a U(4Nf) x U(ANf) 
transformation with U = p\, V = 1. Then the new term in the action takes the form trM. 
In any case, it is clear that AS'efj breaks the U(4Nf) x U(4:Nf) symmetry to U(4Nf). No 
chiral symmetry is left. 

Another way to see the destruction of chiral symmetry is to treat A5 e ff as a perturbation 
on the symmetric action S^°. As we have seen, the latter possesses (4iVj) 2 Goldstone 
bosons; the new term renders all these bosons massive. 

These results show the same pattern as the symmetry breaking in a theory of naive 
fermions due to a fermion mass term. The chiral 17(4) x U (4) symmetry of massless naive 
fermions is broken by a mass term mfipip to U(4). The latter is a vector symmetry, and one 
expects no Goldstone bosons. This is our central result: Domain wall fermions in strong 
coupling behave as naive, doubled fermions with a mass term. 

The same symmetry-breaking pattern emerges from addition of a Shamir mass term || 
to the domain- wall action. This couples the fields on the boundaries s = 0, L 5 according to 
the action 



12 



S s = -pm s ^ *PxL 5 + V'aO (~2^) 



(3.24) 



In the 05 — > limit, this adds a term /15 to the single-site Hamiltonian h. Its matrix element 
between wave functions ip(s) and ip'(s) is 

<V| W> = ^(0)^W + (3.25) 
and it contributes a term to the effective action, 

Sf ff (M) = -N c J2trh s M(x). (3.26) 

Perturbing within the ms = ground states, we write 

tvh s M(x) = tvh s M. (3.27) 
and allow M to rotate within the zero-mode sector as above. As shown in 161, within the 



space of zero modes the operator hs takes the form ms(m — 4)pi, as found above for the 
term induced at strong coupling. 

C. PV bosons restored 



Returning to the full effective action (|2.36j) , we can step quickly through our analysis 



above of the pure fermionic theory. We begin with the case r = 0. The site-site cou- 
pling term in S e s is then symmetric under the graded unitary group U{4NfL^\4NfL^) x 
U(4:NfL 5 \4:NfL 5 ). The h term breaks this symmetry to a vector subgroup, except within 
the space of the fermionic zero modes where a U(4Nf) x U(4Nf) chiral symmetry survives as 
before. The pseudofermion modes are all massive and thus do not add any chiral symmetry 
to S e s- We eliminate the a matrices from the action via spin diagonalization, Eq. ( |3.2| ). 
We assume a saddle point at which the homogeneous mean field commutes with h, 

M{x) = M = Y,\ e ) M M- (3-28) 

e 

Just as h has two blocks, one for the fermionic modes and one for the bosons, so does M e . 
Note that this diagonal supermatrix has no non-commutative elements that are nonzero. 
Since the mean field equations stemming from Eq. ( |3.9| ) decouple the modes of h from each 
other, the solution for M e is unchanged from Eq. ( p.lOj ) ,p| In particular, the Pauli-Villars 
modes have no effect on the fermionic modes. The chiral U(4Nf) x U(4Nf) symmetry of the 
fermionic zero modes is broken spontaneously to U(4Nf), bringing about the appearance of 
16iV| Goldstone bosons in the adjoint representation of U(4Nf). 

Restoring the r-dependent terms in S e g, we perturb the saddle point. The r-dependent 
terms in the effective action are [cf. Eq. ( |3.19| )] 



6 There are two blocks in M £ , but both are given by the formula (3.10). 
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AS eS {M) = — I [-r(-l) :Cj Str7 i {M(x),M(x+j)} + r 2 StrM(x)7 j M(x + j)7 j 



. 3=1,2,3 



+ ^ [-rs 3 (x)Str7 4 [M(x),M(x + 4)] -r 2 StrM(x)7 4 M(x + 4)7 4 ]| . (3.29) 

Again we set M(x) = M and allow M to rotate in the space of zero modes of h. We 
evaluate the supertraces in Eq. (|3.29| ) in the basis of fermion states |0 ± cr),|ecr)ir and boson 
states \eka)s- We use Eq. (|A16|) again in the evaluation of Eq. ( |3.29| ), but now Eqs. ( |A28| ) 
and (|A29|) introduce new, nonzero matrix elements. Dependence on M survives in both 
0(r 2 ) terms in Eq. ( p.29| ), which yield 



AS cS (M) = ^r 2 £ Str I £ 7iM 7i M - 74 M 74 M 



2 

N. 



c 2 ' 

— r 
2 



T ( £< "*l7iM7i|0V) - (OV|7 4 M 74 |oV) 

x(0"V|M|0V). (3.30) 
We evaluate the matrix elements in Appendix to obtain the final result 



jJ2 



AS eS (M) = const. - N c - (J - Ay) ]T tr Pl M. (3.31) 



Only the coefficient has changed from Eq. ( |3.23|) . One might have expected that the well- 



known cancellation between the fermions and the Pauli-Villars bosons affects the low-energy 
sector, but this does not occur. Nonetheless, the heavy superpartners of the mesonic bound 
states must play an essential role in canceling the infinite number of heavy modes, thus 
allowing the above truncation to a low-energy effective theory. 

IV. DISCUSSION 

In this paper we have argued that, at strong coupling, non-Abelian gauge theory with 
domain wall fermions exhibits explicit chiral symmetry breaking but with an enlarged spec- 
trum of pseudo-Goldstone modes characteristic of the doubling phenomenon. We derived 
the effective action of domain wall fermions at the leading order of the strong coupling and 
large- N c expansions in the limit that the lattice spacing of the fifth dimension goes to zero 
and its extent to infinity; therefore our conclusions apply also to overlap fermions. 

Our results also depend on perturbing to lowest order in the Wilson r parameter in the 
neighborhood of a stationary point with maximal symmetry consistent with r = 0. Another 
phase, with different symmetry, might be encountered at finite r. We consider it unlikely 
that moving to a saddle point with lower symmetry will restore any symmetry broken at 
this most symmetric saddle point. Thus it would be surprising if chiral symmetry is restored 
at the true r ^ saddle point. On the other hand, the U(ANf) symmetry characteristic 
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of exact doubling could end up broken; likewise, corrections of higher order in 1/g 2 might 
break this symmetry. We do not know of another signal of doubling that could then be 
sought in the meson spectrum. 

Our analytic results provide evidence for a phase transition from a weak coupling phase 
with good chiral properties to a strong coupling phase where the desired chiral properties 
are lost. This is based on the observation that to lowest order in r there is an explicit 
breaking term identical in form to the quark mass term introduced by Shamir. Thus the 
true Goldstone modes are lost. In the Monte Carlo simulations of domain wall fermions, a 
pattern of increasing violation of chiral symmetry has been observed as the coupling is made 
stronger. The violations of chiral symmetry are suppressed exponentially in L 5 and therefore 
increasing the coupling requires increasing L 5 and consequently the computational burden. 
Our calculation, however, has been done in the L5 — > 00 limit so we should not be seeing 
this problem. In this limit, domain wall fermions are equivalent to overlap fermions. The 
locality of the (effective) 4d Dirac operator is only guaranteed [^3J at moderate couplings 



and therefore the bad chiral properties at strong coupling might be connected with the fact 
the theory is no longer local. 

Since the induced breaking of chiral symmetry is identical in form to that caused by a 
Shamir mass term, it is tempting to consider tuning the latter to cancel the former. This has 



been considered, in the context of finite-L 5 effects, by the authors of ||25|| . One is reminded of 
Wilson fermions, where the pions may be made massless by tuning the hopping parameter. 
Wilson fermions do not regain chiral symmetry by this procedure, except in the continuum 
limit; the masslessness of the pion stems from the proximity of a phase where parity and 
flavor symmetry are spontaneously broken WM. If the analogy is exact, then tuning ms to 



make the Goldstone bosons massless will not restore chiral symmetry in the lattice domain- 
wall theory. Neither would it restore locality to the effective overlap operator. 

Recently, overlap fermions have been studied with a hopping parameter expansion, valid 
for large values of m [E7|. Our coefficients J and y vanish in the limit m — > 00, so that 



chiral symmetry is apparently restored and our results are consistent with those of |27 
Golterman and Shamir [p8fl , however, have argued that a chiral U(Nf) x U(Nf) symmetry 
remains rigorously unbroken for sufficiently large (not necessarily infinite) values of m. Thus 
there should be at least N 2 Goldstone bosons, including a U(l) pseudoscalar. Our results, 
while in agreement with the earlier Hamiltonian study fl6 |, stand in contradiction with those 



of [^8| in the common region of applicability, namely, large g and m. We can only conjecture 
some subtlety in taking the double limit g 2 ,m — > 00. 

Our strong coupling results rest on the analysis of a chiral U (4N fL 5 \4N fL 5 ) super matrix 
model on a 4d lattice. An external field is provided by the Hamiltonian in the fifth dimension. 
Such matrix models are very interesting but notoriously hard to analyze. This supermatrix 
model reformulation of domain wall fermions is not peculiar to strong coupling. At any 
coupling, after integrating out the gauge fields, U(N C ) Yang-Mills theory with domain wall 
fermions must take the form of an effective theory, 



VM(x) exp 



S A {M) -N C Y, Str log M(x) - N c Str hM( 



x 



(4.1) 



expressed in terms of meson-like superfields. This is the ineluctible consequence of local 
4d gauge invariance. The effective theory contains all the mesonic bound states — (4NfL 5 ) 2 
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spin-flavors of quark-antiquark pairs, a like number of pseudoquark-antipseudoquark pairs, 
and their superpartners. The physically interesting limit is 4NfL 5 — > oo, an infinite number 
of "flavors." 

Changing from U(N C ) to SU(N C ) will introduce baryon-like superfields but these are 
most likely irrelevant to the low energy chiral phase structure. With this caveat, the only 
difference between strong coupling and general coupling is that now the 4d action, S^M), is 
no longer restricted to nearest-neighbor couplings, becoming more and more non-local as we 
approach weak coupling. We expect this exact representation for U(N C ) domain wall lattice 
QCD to be in the correct chiral phase at weak enough coupling. The central question of this 
paper, whether the correct chiral symmetry pattern persists at strong coupling, is a question 
concerning the role of the non-locality for a U(n\n) 4d supermatrix model in an external 
field at large n. Our present results give some insight for the nearest-neighbor matrix model 
derived in strong coupling near r = 0. More general analysis of these matrix models can 
help to understand domain wall fermions in the confined phase. 
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APPENDIX A: 



We use the following Dirac matrices for the Euclidean theory: 



Ik = P20~k, k = 1, 2, 3 74 = Pi 

«fc = 747fc = ipz<?k a 4 = l (Al) 

75 = P3 <V = iliflv 



The single site fermion Hamiltonian is the same as that studied in 16 



d 

h = 7475^ + P74, (A2) 

where p = m — 4. We impose the boundary conditions 

(l + 78 )u(L) = (1-tbMO) = 0. (A3) 

on the eigenf unctions of h. We have taken a continuum limit as — > 0, L5 — > 00, with 
L = a^L 5 held fixed. 

Since 7475 = —ip2, there is no spin dependence in h and thus its eigenfunctions take the 
form 

u = ( tyl* ) (A4) 
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where x is an y 2-spinor. The spectrum contains discrete surface states with wave functions 

/ = ±A sinh k(s -L) = ±/oO) (A5) 
g = A smhns = g (s), (A6) 

where k satisfies the eigenvalue condition 

k = /itanh kL. (A7) 

Their energies are 

cosh/tL ^ 

As L — > oo, we have k — > /i and eo — > 0. Taking spin into account, there are four zero-energy 
states for each flavor, which we denote by \0 v a), where rj = ± is the sign of the energy before 
taking L — > oo, and a is a spin label. 
There are also continuum modes, 

f = ±A kS mk(s-L) = ±f e (s) (A9) 
g = A k smks = g e (s), (A10) 

with the quantization condition 

k = fitsmkL (All) 

and energies 

e = ± - 
cos kL 



±Jk 2 + /j 2 . (A12) 



There is a gap to this pseudo-continuum, |e| > fi = m — 4. When 3> /i the solutions of 
Eq. ( [A 11|) approach (2n + 1)tt/2L. We label these states |ecr). Note that we take k > and 
that there are two energies e for each value of k. 
The relation {75, h} = implies that 

75 |0V) = |(TV) and 7 5 |e<r) = | - e, a). (A13) 

Together with orthogonality of the eigenf unctions, this relation can also be used to prove 

dsf f e = / dsg g £ = 0. (A14) 



Moreover, the overlap integral 



dsf (s)g (s) (A15) 

goes to zero exponentially as L — > 00. Using the explicit forms of the 7 matrices, we find 

(0V| 7m |(FV) = (0V| 7fe |«7') = 0, (A16) 
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for n = 1, . . . , 4 and k = 1,2, 3. On the other hand, 



<0V| 74K >= j^M ^ = ' (A17) 



where 



w = £*tM*u*-&wtw < A18 > 

for large L. 

We use these matrix elements to calculate the coefficient matrix in the symmetry- 
breaking term in the effective action (|3.22| ). This matrix is 

(OVI74M74I0V) = ^(0V|7 4 |er)M e (er|7 4 |0' 7 V / ). (A19) 

er 

According to Eq. (|A17 ), the matrix elements of 74 are zero unless rj = sgne and rf = sgne; 
moreover, there is no spin structure in 74 or in M e and so a = r = a'. Thus the matrix is 
diagonal, and we need calculate only (OV I74M74IOV). Using Eq. ( |A13| ), we have 

(O-VI74M74IO-V) = (OVI7475M7574IOV) 

= E(° V l747 5 |ecr)M e (ea|7574|0V) 

€ 

= E( 0V N - e, <r)M e (-e, a^a) 

e 

= E(° V l74|e, ( x)M- e (e,a|74|0V) 

e 

= E(0V|74|e,a>(-M e )(e,a| 7 4|0V>, 

by virtue of Eq. ( |3. 13|) . Thus 

(O-VI74M74IO-V) = -(OVI74M74IOV), (A20) 
and the matrix takes the form J p%. Explicitly, 

J = (0 + +|7 4 M74|0 + +) 

= / rfep(e)(0 + +|7 4 |e+)M e (e+|7 4 |0 + +) 
Jo 

2/x r°° k 2 6y/fc 2 + + 8y/fc 2 + /i 2 + 7 

~ 7T JO (P + p 2 ) 2 k 2 + /i 2 + 16 { ' 

when d — 4. Here p(e) = (L/7r)(e/fc) is the density of states. 

Let us turn to the Pauli-Villars bosons. Their wave functions satisfy Eq. ( |A2fl with 
antiperiodic boundary conditions, 

u(L) = -u(0). (A22) 
The spectrum contains only continuum modes, 
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/ = ±Ae iks = ±f k (s) (A23) 
g = Ae lks = g k (s), (A24) 
where the momenta are quantized as k = ±(2n + l)n/L. The energies of these modes are 



e = i^F + zi 2 , (A25) 

and therefore there is a gap \e\ > p = m — 4. k runs from — oo to oo, and there are two 
energies e for each value of k. The Pauli-Villars fields exhibit no zero mode. 

We will need to know the overlap integrals between the continuum Pauli-Villars modes 
and the fermion zero modes. For large L, 

r L . , , . , . nr 1 



I dsf (s)g k (s)^- ] /^ WT ^(ii + ik) (A26) 
£ ds f k (s)g (s) ~ ^^J—^^ + ih). (A27) 



By using these integrals and the explicit form of the 7 matrices we find[] 

(OVlTfleforOfl = io aa , £^ fc2 ^ 2 [/i (77 - sgn e) + ifc (77 + sgn e)] , (A28) 

(0V| 74 |eA:a% = ^V^ fe 2 + 2 ^ ^ ~ s S ne ) + ^ ("^ + s S ne )] ■ ( A29 ) 
Proceeding to the effective action (|3.30|) , let us first evaluate the matrix 

(OV^-A^oV) = ^(OV| 7j |eA;r) B M eB (eA;r| 7j |0 ,? 'a / ). (A30) 

ekr 

Using Eqs. ( [3.13] ) and (|A28|) one can easily prove that the matrix ( |A30|) takes the form yp%. 



The coefficient is 

y = (o + +| 7 ,m 7j |o + +) 

/oo 
dep{e)\ J2(0 + +h\^+) B M eB (ek+\ ll \0 + +) 
"OO 7„ 



k = ± 



2p r°° k 2 -p 2 6 VP + yu 2 + 8_VP + /i 2 + 7 

~ 7T 7o (F+/i 2 ) 2 F+/i 2 + 16 ' 1 j 

It is very similar to J , and it is independent of the index j. 
We can now evaluate the matrix 

(OVI74M74I0V) = ^(0V|7 4 |eT) F M eF (eT|7 4 |0V) 

er 

+ Y J {Wka\ lA \eT) B M eB {ekr\ lA \W'o'). (A32) 

ekr 

The first sum in Eq. ( |A32|) is the same as that in Eq. ( |A19|) and is equal to J ~p 3 . Using 
Eqs. ( |3.13| ) and ( |A29| ) we find that the second sum equals — 3^3. Changing basis in the 



space of zero modes as in Eq. ( |3.23| ) converts ps — > p\. 



7 We mark the bosonic continuum states with B subscripts and include k in the label because of 
the k — > —k degeneracy; an F subscript is added to fermion states where needed to avoid confusion. 
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APPENDIX B: 



In this appendix we derive the basic identity for the generating function that allows the 
fermionic Grassmann integrals to be replaced by integrals over a unitary matrix field M(x) 
[see Eqs. (|2.17| ) and (p.19 )]. In particular we present the generalization of this well-known 
result to supermatrices |26[ , which we require for the inclusion of the Pauli- Villars fields [see 
Eqs. (EOl and fl2~35D l. 

In the absence of Pauli- Villars fields, we seek the "replacement" of the fermion bilinears 
with matrices M(x) G U(4N f L 5 ), 



1 N ° 

— Y ip aa ip af3 



" c a=l 



The essential identity allows the replacement of the generating function 



Z(J) = J dip dip exp (-JVctr JM) = (det J) 



by the U(n) group integral 

I(J) = 



dM 



(det M) 



N c 



exp (-JVctr JM) , 



(Bl) 



(B2) 



(B3) 



up to a multiplicative constant, where dM is the Haar measure for U(n) and n = 4JV^L 5 . 
This identity was proven by Kawamoto and Smit |TU[ and it follows immediately from a 



similar identity for the one-link integral given earlier by Creutz [19 



The generating function that arises in the theory with Pauli- Villars fields is 



Z(J) = dip dip d<p d(P exp (-JV C Str JM) = (Sdet J) 



N c 



(B4) 



where M is the supermatrix given in Eq. ( |2.26| ). The generalization to supermatrices of 
Eq. O is 



7(J) 



dM 



(Sdet M) 



exp(-JV c Str JM) 



(B5) 



where dM is the invariant Haar measure for the supergroup U{n\n). The latter is the 
group that leaves invariant the norm J2i z i z i + Si^i^i m complex superspace, where 6i are 
Grassmann variables and Zi are c-numbers. 

We choose definitions of the supertrace and superdeterminant that differ from the con- 
ventional choice by Str — > —Str and Sdet — ► 1/Sdet in order to simplify the comparison 
between expressions before and after the including the Pauli- Villars fields.0 To be explicit, 
we write 



8 See Efetov 
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with the supervector index in the order (9 1 , 9 2 ■ ■ ■ 9 n , z±, z 2 • • • z n ), i.e., Grassmann variables 
first. Then 

StrM = trA-tr.D, (B7) 

and 

The proof of our identity relies on the following observations. Since the integral measure 
is compact and the integrand is free of singularities, the integral I (J) is holomorphic in the 
components of J. In particular, if we write 

j = | J ^ J ^ | (B9) 

J<fnji J<f><f> 

in block form and replace — > zJ^ and — > z -1 J^, we can expand I (J, z) in a Laurent 
series, 

oo 

I(J,Z)= CmZ™, (BIO) 

m=— oo 

where I(J) = I (J, z — 1). There can be no dependence on z*. 

Next we factor the Haar measure for U(n\n) = U(l) x C/(n|n)/C/(l) as = dcpdM. 
The explicit parameterization for this £7(1) factor is 

M= ( C e-^J- (BU) 
where M = M{<j> = 0) is special unitary, Sdet M = 1. Note that Sdet M = e*™^ Sdet M = 
Introducing the source rotated under U(l), 



J = 



ze ^ 

7 - 7 -l„-!^/2 T 

y z e i 



(B12) 



the group integral becomes 

I(J,z) = J d<pdM -J^exp (-7V c Str JM) = J d<f)e- iN ^F{J). (B13) 

Moreover, invariance of the integral dM over SU (n\n) = U (n\n)/U(l) implies that the kernel 
F(J) is invariant under right and left rotations, 
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F(J) = F{UJV ] ), 



(B14) 



with U, V G SU (n\n). Thus F(J) must be a function of the invariants Sdet J and Str (JJ^) 1 . 
Holomorphy rules out the trace terms Str (JJ^) 1 because of their explicit dependence on z*. 
Then using Sdet J = (ze^) n Sdet J, we see that the integral over projects out a single 
term in the Laurent expansion, c m for m = N c n. Hence 



I(J, z = l) = c m (Sdet J) 



(B15) 



proving the identity. 

In the remainder of this appendix we shall prove Eqs. (|2.13|) and ( 2.311) , which are of 
paramount importance in the derivation of the effective actions ( |2.22| ) and ( |2.3(j| ). The proof 
of Eq. ( |2. 13|) was given in [10| and we shall briefly sketch it here in order to extend it to 
the case of Eq. ( 2.31 ), in which one has also to deal with the Pauli-Villars fields. Using 
Eq. ( |2.9| ), we express the action ( |2.13| ) as a sum over one-link integrals, 



S l {il>,if>) = Y,w(A{x),A(x)) =J] log/ dU expTr \A„(x)U + U ] A 



X,fl 



(B16) 



where A and A were defined in Eqs. ( |2. 10|) and ( 2.11 ). The one-link integral appearing in 



Eq. (|B16|) was computed in the large- iV c limit in pl|j22|, with the result 



w 



(A, A) = N 2 C I I - c + 2N C J2(c + x a ] 



1/2 



2iV 2 

c a, i 



^ log {c + xaf^ + ic + x,) 1 ' 2 ). (B17) 



1/2 



Here the eigenvalues of AA/N 2 and c is defined implicitly by 

1 



2N, 



-1/2 



c a 



By expanding in x a , Eqs. ( |B17| ) and ( |B18| ) can be rewritten as a series in 
By rearranging the fermion fields in A and A we obtain 



Tr 



AA 

N*. 



-tr 



-lA(x) 



(B18) 



(B19) 



(B20) 



where X(x) is defined in Eq. (|2.12| ). Thus for any function that admits a series expansion 
we can write 



^E/(*-) = /(o)-F tr 



/(-Ja)-/(o) 

and using Eqs. ([5T7D -([B2lD we obtain Eq. (gl3|) . 



(B21) 
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By the same method we can prove Eq. (|2.31|) in the presence of the Pauli-Villars fields. 
The proof relies on the fact that Eqs. ( B17|) -( B21 ) still hold but the trace tr should be 
replaced by the supertrace Str. Moreover A and A are now defined by Eqs. ( |2.29| ) and ( j2.30| ), 
and M. by Eq. Q2.26Q . One can easily verify that in presence of Pauli-Villars fields one has, 
instead of Eq. flB~2U 




-Str 



(B22) 



Equation ( |B22j ) allows the expansion of Eq. (|B17|) in powers of x a in the supermatrix case 
as well. One can thus prove that the action can be completely expressed in terms of A as in 
Eq. ([2.31 ) , where now A is a bilinear function of the supermatrix Ai given in Eq. ( |2.26| ) . 
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FIGURES 
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FIG. 1. The function F(X). The dotted line is F = A/4. 
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FIG. 3. The integrals J (dashed curve) and y (dotted curve), and the combination J — \y 
(solid curve), as functions of fx = m — 4. 
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